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ABSTRACT 

The analysis of the four-colour maps of galaxies in the Hubble Deep Field has revealed 
in the sample of 0.4< z <1 early- type field galaxies the existence of ellipticals with a 
predominantly old coeval stellar population. However, there is another, unexpected, 
category of HDF early- type galaxies, in which the galaxy core is significantly bluer 
than the outer regions. We demonstrate that these colour gradients are predicted by 
the multi-zone chemodynamical model for the evolution of elliptical galaxies. 

We suggest that the colour gradient could be used as a chronometer for the evo- 
lution of elliptical galaxies: galaxies younger than a few Gyr exhibit cores bluer than 
the surrounding galaxy, due to ongoing star formation while more evolved galaxies 
have redder cores, due to metallicity gradients increasing toward the centre. 

Key words: cosmology: observations - galaxies: elliptical - galaxies: evolution - 
galaxies: formation- galaxies: ISM - galaxies: starburst 
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1 INTRODUCTION 

A surprising result of the analysis of the four-colour 
((730oi345oV606-^8i4) niaps of the Hubble Deep Field (Abra- 
ham et al. 1999, hereafter A99; Menanteau, Abraham & Ellis 
2000) is the existence of early-type field galaxies, with posi- 
tive colour gradients in the centre (galaxy nucleus bluer than 
the outer surrounding region) evidencing recent nuclear star 
formation in early-type galaxies with 0.4< z <1. 

Although these bluer core spheroids (BCSs) could be 
explained within the scenario of structure formation via hi- 
erarchical merging, in which the recent nuclear star forma- 
tion would be the result of mergers leading to the formation 
of an elliptical galaxy, the high degree of symmetry in the 
images requires, additionally, that, after the merger, the sys- 
tems had had time to relax to a spheroidal geometry. 

On the other hand, the appearance of these systems is 
exactly that predicted by the chemodynamical model for 
evolution of elliptical galaxies (Friaga & Terlevich 1994; 
Friaga & Terlevich 1998, hereafter FT98). The chemody- 
namical model is a single collapse model, diiTering, however, 
from the monolithic collapse model, in that it is a multi- 
zone model, in which there is no coordination in the star 
formation in the several zones of the model. This allows 
that, locally, the time scales for star formation could be of 
the order of 0.1 Gyr in the central regions of the galaxy, 
leading to a super-solar [Mg/Fe] ratio as observed in lumi- 
nous ellipticals, while it takes ~ 1 Gyr for the main body of 
the elliptical galaxy to be formed. An important ingredient 
of the chemodynamical model is the persistence of a central 



cooling flow for 2 — 3 Gyr, feeding star formation in the core, 
which would account for the core being bluer than the rest 
of the galaxy, as in the BCSs found in the HDF. 

Within the frame of the chemodynamical model, the 
BCSs at intermediate redshifts would be one of the mani- 
festations of the general evolution of spheroids. At higher 
redshifts, the chemodynamical model for the evolution of 
spheroids predicts a link of the BCSs with Lyman break 
galaxies (Steidel 1995, 1996) and the scarcity of passively 
evolving elliptical galaxies (Zepf et al. 1997, Barger et al. 
1999). At z 3, the first ^ 1 Gyr of our model galax- 
ies shows striking similarities to the Lyman break galax- 
ies (LBGs): intense star formation, compact morphology, 
the presence of outfiows, and significant metal content. Our 
investigation of the nature of the LBGs using the chemo- 
dynamical model for spheroids ( Friaga & Terlevich 1999, 
hereafter FT99) supports a scenario in which LBGs are the 
progenitors of the present-day (age of 13 Gyr) moderately 
bright (0.051/* — 1.41/*) spheroids. Over the z > 1 range 
in general, the continuing (for ~ 2 — 3 Gyr) star formation 
in the centre of ellipticals predicted by the chemodynamical 
model could explain the relative absence of very red galaxies 
in deep optical and near-infrared surveys, as expected from a 
population of passively evolving elliptical galaxies (Jimenez 
et al. 1999). The BCSs would indicate the persistence of 
the same continuing central star formation at redshift < 1. 
It should be noted, however, that most of the mass of the 
stellar population is formed during the first Gyr; typically 
half of the present stellar mass of the galaxy is formed in 
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~ 0.5 Gyr, and 5% — 30% during the extended lower level 
central star formation phase. 

In this paper, we follow in detail the evolution of colour 
gradients in spheroids within the scenario of the chemody- 
namical model and investigate the relation of the continu- 
ing central star formation activity to the BCS population. 
Section 2 gives a description of the chemodynamical model. 
Section 3 presents the predicted evolution of colour gradi- 
ents in spheroids. We show that, for a wide range of masses, 
there are 3 phases: 1) an earlier phase, with relatively flat 
colour gradients; 2) the bluer core phase; and 3) the redder 
core phase. Section 3 also discuss the possibility of using the 
transition from the bluer core phase to the redder core phase 
as an evolutionary clock in cosmological studies. Finally, in 
section 4, we identify the BCS soon at 0.4 < z < 0.7 as 
sub-L* spheroids formed at 2; > 1 with continuing central 
star formation or recurrent bursts of central star formation. 
Section 4 summarises our conclusions. 

We adopt in this paper Ho ~ 70 km s~^ Mpc~^, in ac- 
cordance with the value obtained by HST Key Project on the 
Extragalactic Distance Scale, Ho = 71 ± 6 km s~^ Mpc""', 
or, including a possible metallicity dependence of the Chep- 
heid period- luminosity relation, i!fo = 68±6 km s~^ Mpc~^ 
(Mould ot al. 2000). We also adopt a flat Q^n + = i 
cosmology. 



2 THE CHEMODYNAMICAL MODEL 

In this work, the evolution of spheroids is modelled with the 
aid of the chemodynamical model of FT98. The chemody- 
namical model combines multi-zone chemical evolution with 
1-D hydrodynamics to follow in detail the evolution and ra- 
dial behaviour of gas and stars during the formation of a 
spheroid. Within each spherical zone, in which the model 
galaxy has been divided for the 1-D liydrodynamical cal- 
culations, the evolution of the abundances of six chemical 
species (He, C, N, O, Mg, Fe) is calculated by solving the ba- 
sic equations of chemical evolution. The star formation and 
the subsequent stellar feedback regulate episodes of wind, 
outflow, and cooling flow. The knowledge of the radial gas 
flows in the galaxy allows us to trace metallicity gradients, 
and, in particular, the formation of a high-metallicity core 
in ellipticals. 

FT98 built a sequence of chemodynamical models re- 
producing the main properties of luminous elliptical galax- 
ies. The calculations begin with a gaseous protogalaxy with 
initial baryonic mass Mg- FT98 investigates models with 
Mg between lO" and 5 x 10^^ M©. In order to study LBGs, 
FT99 have subsequently extended the mass grid of the mod- 
els of FT98 downwards to 5 x 10^ Mq. Intense star for- 
mation during the early stages of the galaxy builds up the 
stellar body of the galaxy, and during the evolution of the 
galaxy, gas and stars exchange mass through star formation 
and stellar gas return. Owing to inflow and galactic wind 
episodes occurring during the galaxy evolution, its present 
stellar mass is ~ 15 — 70% higher than Mq- Gas and stars 
arc embedded in a dark halo of core radius rh and mass 
Mh (we set Mh ~ 3Mg). The models are characterised by 
Mg, Th, and a star formation prescription. The SFR is given 
by a Schmidt law vsf oc p^^'^ (p is the gas density and 
vsF = SFR/p is the specific SFR). Here we consider the 
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Figure 1. Evolution of the colour C/300 — S450 (rest- frame) 
through several annular apertures for the models with Mg = 
2 X 10" Mg (labeled Model 2), with Mg = 2 x IQi^ Mq (Model 
20), with Mg = 10^2 M© (Model 10), and with Mq = 5 x 10^ 
Mq (Model 0.05). 

standard star formation prescription of FT98, in which the 
normalization of v is vo = 10 Gyr~^ (in order to repro- 
duce the super-solar [Mg/Fe] ratio of giant ellipticals) and 
xsF ~ 1/2 (with the exception of one model, for which 
xsF = 1/3). The stars form in a Salpeter IMF from 0.1 
to 100 Mq. a more detailed account of the models can be 
found in FT98. 



3 EVOLUTION OF COLOURS GRADIENTS 

We have used the spectrophotometric models of Bruzual & 
Chariot (1998) to compute the SED for each galaxy zone. 
This is done in a self-consistent way using the metallicity 
given for each zone and each time by the chemodynamical 
model and for a consistent IMF. Then, the SED is inte- 
grated over several apertures, redshifted, reddened by the 
Lyman a forest opacity (Madau 1995), and convoluted with 
the filter transmission curves. Unless otherwise stated, AB 
magnitudes are used throughout this paper. 

It is helpful for the understanding the evolution of the 
colour gradients to consider first the evolution of the colours 
in the rest-frame of the galaxy. Figure 1 shows the evolution 
of the t/300 — B450 colour, more sensitive to on-going star 
formation, as seen through several annular apertures, for 
several models, with Ma ranging from 5 x lO'^ to Mg = 
2 X 10^^ Mq. The evolution of model with Mg = 2 x lO" 
(the fiducial model of FT98, or model 2 —following FT98 
we label the models according Mg in units of 10^^ Mq) is 
typical of models covering a wide range of masses (from less 
than 10^° up to 10^^ Mq). Its evolution exhibits 3 distinct 
epochs regarding the colour gradients: 1) An early phase, 
during which the galaxy as a whole is blue, due to wide 
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spread intense star formation. The region with the highest 
SFR's extends out to r ~ 3 kpc with httle radial gradient 
of specific SFR, which results in a colour gradient flat in 
this region. The very central region can be even redder than 
its immediate surroundings during this phase, due to the 
very fast chemical enrichment for r < 1 kpc (where, the 
metallicity, at 0.1 Gyr, of the stellar population is 0.6 Zq) 
2) After 3 x 10* yr, the SFR decreases from outside in, but it 
persists for more than 1 Gyr in the core (r < 1 kpc), which 
is significantly bluer than the rest of the galaxy: that is the 
bluer core phase, during which the galaxy would be seen as a 
BCS; 3) The last phase {t > 1 Gyr) is the redder core phase, 
during which the level of star formation is low in the core, 
which is redder than the outer parts of the galaxy, due to 
its high metallicity. The dips in the evolution of U^oo — B450 
indicates bursts of star formation. Note that the transition 
from bluer core phase to the redder core phase is rapid, 
which suggest the use of this transition as a time-marker. 
The transition occurs later if we consider the surrounding 
2 < r < 3 kpc region, because in the outermost radii (5 < 
r < 10 kpc), the lower metallicities imply bluer colours. 

The evolution of the model with Mq = 2 x 10^^ Mq 
(model 20 of FT98) illustrates another sequence of phases 
that occurs only for very massive systems {Mq ^ W^^ Mq). 
In this mass range, the galactic wind occurs late, and the 
inner cooling flow stops for a short time after which it reap- 
pears and eventually becomes a massive, global cooling flow 
extending over 100 kpc. The bluer core is clearly present 
only in the early stages of the evolution (till ~ 0.3 Gyr). 
After this time, due to the complex inner cooling flow be- 
haviour, which extends over the whole r < 10 kpc, several 
uncoordinated star formation episodes rend the colour gra- 
dients ill-defined. For t 1.5 Gyr, the core is redder than 
its surroundings, although the transition of a bluer core to 
a redder core is not smooth, and the evolution of the colour 
of the core is subject to some blueing at ~ 3 Gyr due to 
a weak, late, central star formation episode. Finally, in late 
times t ^ 10 Gyr, the establishment of a massive, global 
cooling flow recovers a bluer core. 

Note that, although the model with Mq = 10^^ M© is 
also massive enough to develop a cooling flow in t > 10 Gyr, 
(visible as a dip in the U-Mm — -B450 profiles after this time), 
during its first 10 Gyr it exhibits the three-phase evolution of 
the colour gradients, with the characteristic cleax separation 
between the bluer core and the redder core stage at ~ 1 Gyr. 

3.1 A typical ~ L* galaxy 

The fiducial model of FT98 [Ma = 2xl0" M© and r/. = 3.5 
kpc ) has a present-day stellar mass, 2.4 x 10^^ M©, corre- 
sponding to Lb ~ 0.7L* , being typical for a present-day 
moderately luminous elliptical galaxy. Figure 2 shows the 
evolution of the colour gradient for this model. The quan- 
tity AC is defined as the colour at an annulus 5 < r < 10 
kpc minus the colour in the central kpc (r < 1 kpc). Pos- 
itive values of AC therefore correspond to positive central 
colour gradients, and indicate a core bluer than the sur- 
rounding galaxy. The 5 < r < 10 kpc dimension chosen for 
the annulus representing the non-nuclear part of the galaxy 
falls inside a typical effective radius of a luminous elliptical 
(~ 10 kpc for a L* elliptical). This region is easily resolv- 
able by large ground-based telescopes, since 1 arcsec corre- 
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Figure 2. The evolution of the colour gradients AC of the 
colours: Veoe — ^814 (heavy solid lines), B450 — Veoe (light solid 
lines), and C/300 — -B4B0 (dotted lines), for the model with 2 x 10^^ 
MQand for several cosmologies and two epochs (z^gf) of galaxy 
formation. 

sponds to 8.0 (9.1, 7.4) kpc for ifo = 70 km s"^ Mpc"^ 
and Q.m = 0.3 (= 0.12, = 0.45) (assuming a flat cosmol- 
ogy flm + f^A = 1)- At the same time, the r < 10 kpc re- 
gion avoids the oc (1 + z)^* Tolman dimming of the surface 
brightness which makes it difficult to detect the outermost 
regions of the galaxy, if it is at high redshifts. As we can see 
from Figure 2, the changes in the central colour gradient are 
more dramatic for an epoch of galaxy formation zgf = 5, 
but they are very noticeable for zgf = 1- If the gala:xy is 
formed at a recent epoch {zgf = 1), the bluer core is more 
noticeable in the Veoe — Isi4 colour, while for higher redshifts 
of galaxy formation , the C/300 — -B450 colour would trace bet- 
ter the continuing star formation in the galaxy centre. For 
instance, for zgf = 5, the galgixy would have an bluer core 
for 2; < 3 (when it is in the phase of central continuing star 
formation) only in the C/300 — B450 colour (with a maximum 
around 2 = 2 for zgf = 5), while for the B450 — 14o6 and 
V6O6 — ^814 colours, the core is redder than the surrounding 
galaxy. 

The bluer core phase may be characterised by Zrev, the 
redshift at which there is the reversal from blue core (due 
to central star formation) to red core (due to metallicity 
gradient). The values of Zrev are exhibited in Table 1 for the 
models, colours and cosmologies used in this paper. 

The colour Veoe — ^814 (in which we can follow more 
closely the evolution of the colour gradient) has Zrev = 0.764 
for a ftm = 0.3, = 0.7 cosmology and a recent galaxy 
formation epoch {zgf = 1)- The maximum of Veoe — .^814 
is AC = 0.414 at z = 0.859. There is also a secondary 
maximum of the central star formation at z = 0.799 (AC = 
0.285). At the maximum of star formation, this galaxy would 
be very bright (/814 = 17.74) and it would be easily de- 
tectable in the HDF. Since no BCS as bright as this has 
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Figure 3. The same as Figure 2, but for the model with 2 x 10^^ 
Mq. a Clrn = 0.3, Qa = 0.7 cosmology is adopted here. 



Figure 4. The same as Figure 3, but for the model with 2 x 10^^ 
Mq and xsf = 1/3. 



been detected in A99 sample, we can tentatively conclude 
that the typical BCSs are not ~ L*-ellipticals formed at 
ZGF ~ 1. Turning to higher zgf's {zgf ~ 5), the blue core 
is more pronounced in the U300 — B450 , in which the max- 
imum colour gradient occurs at z = 2.012 (AC = 1.851). 
Note, however, that the model galaxy would be practically 
invisible in the C/300 band {U300 = 29.51 for the whole 
galaxy), while still detectable in other bands (-B450 = 25.66, 
Veoe = 24.55, and 78i4 = 23.71). As a matter of fact, it 
would be classified as a LBG, because its colours B460 — Veoe 
=1.11 and Z/soo — B460 =3.85 satisfy the criteria for se- 
lecting LBGs at 2 < z < 3.5: -B450 — Veoe < 1.2 and 
t/300 - 5450 >B450 - Veoe +1.0 (Dickinson 1998). Note that 
for a Qrn = 0.12, = 0.88 cosmology, the luminosity dis- 
tance would be larger and the galaxy fainter: U300 = 31.05, 
B450 = 27.97, V6O6 = 26.23, and hu = 24.99 at « = 2. Then, 
B450 — Veoe =1.74 and it would no longer be identified as a 
LBG. 



3.2 A massive central galaxy in a cooling flow 

Figure 3 shows the evolution of the colour gradient for the 
FT98 model with Mg = 2 x 10^^ Mq, which has a present- 
day stellar mass of 2.4 X 10^^ Mq. Note that the late cooling 
flow episode in this model would require enough gas supply 
through the outer boundary of the galaxy, a condition that 
would bo mot if the galaxy would be located at the centre of 
a relaxed cluster of galaxies. As a matter of fact, the model 
with Mg — 2 X 10^^ Mq could represent an elliptical at the 
bottom of the potential well of a cluster of galaxies collecting 
the gas of a cooling flow formed in the intra^cluster medium, 
since the most massive central galaxies of cooling flows have 
masses of this order. 

On the other hand, the early cooling flow episode could 
be maintained by the gas within the dark mass halo of the 
protogalaxy itself and thus would not require that the galaxy 
be located in a cluster, or, more properly speaking, that the 
protocluster containing the galaxy be well-developed. 

In this model, the first episode of star formation is over 
a more extended region so the central colour gradient is 
smaller. In addition, for a galaxy formed at zgf ~ 5, the 
reversal blue core/red core is followed by a contra-revorsal 
red core/blue core due to later central star formation, which 
would correspond to the present-day central star formation 



activity seen in a few central dominant galaxies in cooling 
flows (McNamara 1997). 

At late times, Model 20 shows starbursts in its core fed 
by the late cooling flow. As a consequence, as we approach 
« ~ 0, AC varies wildly, being positive most of the time 
(star formation in the core prevails). 



3.3 Recurrent star formation in ellipticals 

One particularly interesting model is the model 2(1/3) of 
FT98 because it exhibits recurrence of late central inflow 
episodes. This model has Ma = 2 x lO'^'^ Mq, as the fiducial 
model, but xsf = 1/3 instead of xsf = 1/2 — xsf is 
the index of Schmidt law describing the dependence of the 
specific star formation rate on the gas density vsf <x. p^^'' 
(see end of Section 2). As the SN I rate decreases, the wind 
stalls, the gas starts to accumulate in the galaxy core, and 
the radiative losses drive a new central cooling flow. The late 
central inflow episodes are brief and quickly put out by the 
resulting burst of star formation. The duration of the late 
cooling flow episodes is a few times lO'^ yr, and, since the 
separation between them is typically a few Gyr, their duty 
cycle is ~ 1/100. The turning on and off of the central inflow 
could corrospoiid to the "active" and "inactive" states of the 
core in both the super-massive black hole model of AGN 
and the starburst model of AGN. In addition, the length 
of the duty cycle of the central inflow increases with time 
and the amount of material deposited decreases, implying a 
steep decrease of the time-averaged luminosity of any central 
activity powered by the infalling material. For the model 
2(1/3), the second central inflow occurs at t = 4 Gyr and 
involves 1.5 x lO'' Mq, and the third one (8.4 x lO*' Mq) 
happens a,t t — 9.5 Gyr. Each late central inflow episode is 
accompanied by enhanced star formation inside the central 
kpc, the signature of which is a temporary blueing of the 
galaxy core. 

These late episodes of star forniatiou of the model 
2(1/3) imply that the blue core persists even at redshifts 
lower than in the fiducial model. This can be seen from 
the lower Zrev of the model 2(1/3) in comparison to that of 
model 2. For zgf ~ 1, the first central cooling flow episode 
is more conspicuous in the Veoe — ^814 , as was the case for 
model 2. The maximum of the Veoe — .^814 colour gradient 
occurs at 2 = 0.853 when AC = 0.473 (in the B4B0 — Veoe 
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colour, AC — 0.222). However, the second central cool- 
ing flow episode is more noticeable in ^450 — Veoe , cor- 
responding to the local maximum AC = 0.179 at z = 0.744 
(AC = 0.114 in Veoe — ^8i4 ). The third inner cooling flow 
episode has a more pronounced signature in t/300 — B450 
, but the associated starburst is so weak that it does not 
change the redder core to a bluer core, so AC = —0.012 at 
2 = 0.681. The metallicity gradient effect dominates over 
the impact of the star formation, as we can see from nearly 
constancy of the Veoe — -^8i4 colour gradient (AC = —0.071 
at z = 0.681). 

3.4 Colour gradients as cosmological probes 

Star formation activity in high redshifts galEixies has already 
been proposed to be used to determine the cosmological con- 
stant (Qa) and the matter mass density of the Universes 
(Qm) (Melnick, Terlevich & Terlevich 2000). We suggest 
here the use of colour gradients as a cosmological probe. The 
time-scales for the evolution of the central colour gradient in 
spheroids can potentially provide an independent clock for 
measuring the age of the Universe. Figures 2 also exhibits the 
results for three different cosmologies, illustrating the sen- 
sitivity of the evolution of the colour gradient to the value 
of fim. We consider a flat cosmology flm + = 1 with 
0.12 < flrn < 0.45, as favoured by the simultaneous con- 
straints of the first acoustic peak of the angular power spec- 
trum of the cosmic radiation background measured by the 
BOOMERang (de Bernardis et al. 2000) and the redshift- 
magnitude relation derived from high redshift supernovae 
surveys (Perlmutter et al. 1999, Riess et al. 1998, Schmidt 
et al. 1998) . These latter works have already found evidence 
for an accelerating Universe in which matter would provide 
only ~ 30 % of the critical density for a flat Universe, with 
the rest being accounted for by the cosmological constant 
A. Therefore, the three cosmologies adopted here consider 
Qm = 0.3, = 0.12 and = 0.45. 

An additional constraint is the age of the Universe 
given by each cosmology: 11.98, 13.47 and 17.09 Gyr, for 
a„ = 0.45, = 0.3 and = 0.12, respectively. The « 10% er- 
rors on Ho prevent us from using values of Ho signiflcantly 
lower than 70 km Mpc^^, so the age of the Universe 
implied by f2m, = 0.45 would be dangerously low in compar- 
ison to the ages of the oldest globular clusters, ~ 11.5 — 14 
Gyr, derived in the post-HIPPARCOS age (D'Antona, Caloi 
and Mazzitelli, 1997; Salaris, degl'Innocenti and Weiss, 1997; 
Pont et al. 1998; Chaboyer et al. 1998; Carretta et al. 2000). 

We should also take into account the constraints on the 
age of the Universe from high redshift old galaxies (Lima 
& Alcaniz 2000). One of the strongest lower limits on the 
age of the high redshift Universe comes from the very red 
galaxy 53W069 at z=1.49, with a minimal stellar age of 4.0 
Gyr (Dunlop 1999), which, for Ho > 60 km s"^ Mpc"\ 
implies Clm < 0.5. 

A useful quantity that could be used to set the colour 
clock of the evolving spheroid and distinguish between sev- 
eral cosmologies and investigate the epoch of galaxy forma- 
tion is Zrev The cosmologies with smaller Qrn's (larger ^a's) 
have larger lookback times, and therefore, higher Zrev's. 

Moreover, the models with masses around that of an 
Z/*-elliptical progenitor (Mq = 10^^ — 10^^ Mq) have Zrev 
remarkably similar, which guarantees a fair accuracy for the 
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Figure 5. Top panel: evolution of the colour gradients AC of 
the colours: Veo6 — ^814 (heavy solid lines), S450 — Veoe (light 
solid lines), and C/300 — S450 (dotted lines), for the model with 
Ma = 5 X 10^ Mq, for the flm = 0.3, Ua = 0.7 cosmology and 

an epoch of galaxy formation zqp = 1. Here AC is defined as the 
colour seen inside the annulus 2 < r < 3 kpc minus the colour 
in the central kpc (r < 1 kpc). Lower panel: evolution of the 
^606 — ^814 colour through several annular apertures. The dot- 
dashed line represents the Veoe ^ ^8i4 colour through a 2 < r < 3 
kpc aperture for zqf = 1-2. Also given the evolution of /814 for 
the whole galaxy. 



colour gradient clock. For a Q.„i = 0.3, Q,a = 0.7 cosmology, 

the model with 10^^ M0 has Zrev ~ 0.745 {zqf = 1 and 
Veoe — /814 ) and Zrev — 1.159 {zgf — 5 and Uioo — B450 ). 
The respective values for the model with 10^^ Mq are 0.806 
and 0.970, and, for model 2, 0.764 and 1.263. 

Model 20 is not particularly useful as a colour clock, due 
to the massive late cooling flow, which prevails at the present 
time. It is worth to note that if the model is old enough (zgf 
high and Q,m low) , there is central star formation driven by 
the late cooling flow. Therefore, as we can see from Table 
1, the model galaxy has a bluer core at the present day for 
cosmologies and zgf's that imply larger lookback times till 
the epoch of galaxy formation, U^gf, and in this case Zrev 
is not defined. 



3.5 Bluer Core Spheroids 

It is interesting that the lower mass models of FT99, in- 
tended to describe typical LBGs with velocity dispersions 
around 70 km s~^ have properties that are similar to those 
of typical BCSs of A99 (for instance, BOS HDF-2.264 at 
z=0.478). In particular, the models with Mg = Mq 
and 5 x 10® Mq have prolonged central star formation ex- 
tending for more than 3 Gyr (see Figure 1 of FT99) . Because 
of this, they could be seen as BCSs over a significant fraction 
of their lifetime. These objects, if they have been formed at 
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Table 1. The redshift of reversal from blue core to red core 
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zgf ~ 4 — 5, would be LBGs at z ~ 3, but if zqf ~ 1 
instead, they would be classified as BCSs at z ^ 0.7. 

Figure 5 illustrates the evolution of the colour gradients 
of the FT99 model with 5 x 10^ Mq (hereafter model 0.05), 
assuming zgf = 1- Model 0.05 would be seen as a a = 55 
km s~^ LBG at an age of 0.5 Gyr, and as a present-day (age 
of 13 Gyr) 0.05L* spheroid. In this model, the central cooling 
flow is extinguished only at 3.63 Gyr. By comparison, in 
the fiducial model 2, the inner cooling flow ends at 1.80 
Gyr. In order to compare our results with observations with 
the same resolution as those of HST as in A99, now AC is 
defined as the colour inside the annulus 2 < r < 3 kpc minus 
the colour in the central kpc, since the WFPC-2 pixel size 
of 0.1" corresponds to 0.6 kpc at z = 0.5 for a Q,m = 0.3, 
CIa ~ 0.7 cosmology. The resolution is even better in the 
colours V — I maps of A99, which have a drizzled pixel scale 
of 0.04", allowing, for instance, a clear definition of the r ~ 
0.25" blue core of the BCS HDF-2.264. In addition, in this 
section, we use Vega magnitudes, as in A99. 

As we see from Figure 5, for z > 0.6, our model galaxy 
would be easily detected as a BCS in the Veoe — ^8i4 colour, 
which is the more reliably determined colour in A99 (Menan- 
teau et al. 2000). At z = 0.7, AC=0.472 and Isi4 = 21.52, 
while the colour gradient is fiat in the outer part of the 
galaxy (Veoe — ^814 (2-3 kpc)=1.354 and Veoe - -^814 (5-10 
kpc)=1.311). Note that the galaxy is detectable as a BCS at 
lower redshifts in S450 — Veoe and U300 — S450 [zrev = 0.583, 
0.496 and 0.395 in V^joe - hiA , -B450 - Veoe and Uaot) - -B450 
, respectively). Model 0.05 was not tailored to reproduce 
the BCS HDF-2.264, but the similarities of the properties of 
the model to the object suggest that, if HDF-2.264 is typi- 
cal of BCSs, they would-be z ~ 0.5 counterparts of sub-L* 
spheroids, since /814 = 21.73 of the model at z — 0.5 is 
comparable to /814 = 21.66 (Menanteau et al. 2000) — /814 
of our model varies within the range 21.2 — 21.8 between 
z = 0.4 and z = 0.8. 

The colours of the outer region (2 — 3 kpc) constrain the 
average age of the stellar population of the galaxy. Within 
the scenario of our model the colours of the outer parts 
of HDF-2.264 suggest zgf somewhat higher than 1.0. As a 
matter of fact, model 0.05 has Veoe— -?^8i4 (2-3 kpc) = 1.171 at 
z = 0.5 for Zgf = 1, while it has Veoe — .^8i4 (2-3 kpc)=1.223 
at z = 0.5 for zgf = 1-2, thus reproducing well the red 



outskirts of HDF-2.264, for which the reddest pixels have 
Veoe — .^814 ~ 1-3, with an error of 0.07 mag (see Figure 4 of 
A99). 

On the other hand, the recent central star formation at 
z = 0.478 of HDF-2.264 could be explained by some sec- 
ondary central infall similar to that of model 2(1/3). Note 
that, at t = 1 Gyr, the model 0.05 has a fully developed 
galactic wind with 1200 km s~^ at r = 10 kpc. However, at 
= 3 Gyr, the galactic wind has weakened (wind velocity of 
270 km s~^ at r = 10 kpc) and, in the inner regions there 
is only a highly subsonic outflow (20 km s~^ at r = 1 kpc). 
The gas accumulated in the dark halo of the galaxy brakes 
the outflowing gas, so it is possible that, in some cases, the 
core accumulates enough gas returned by the evolving stars 
to trigger late central star formation episodes, (this happens 
in model 2(1/3) but for a ~ L* galaxy). In this connection, 
we could explain the scarcity of BCSs in galaxy cluster envi- 
ronments (Menanteau ct al. 2000) as due to the stripping of 
the gas in the dark halo of the spheroid, thus preventing very 
prolonged (> 3 Gyr) cooling flow-fed central star formation 
as well as late recurrent central star formation episodes. 



4 CONCLUSIONS 

The early-type galaxies with blue cores of A99 are very sym- 
metrical, showing no morphological peculiarity. A99 points 
out that the ~ 40% fraction they flnd for the early-type 
systems with recent star formation is consistent with the 
predictions from hierarchical models. However, in the hier- 
archical model scenario, after a merger event, the resulting 
system has to relax quickly enough in order to be classified 
as an early-type galaxy. 

A alternative scenario to the merger picture would be 
that proposed by the multi-zone single-collapse model, in 
which late gas central inflow would feed star formation in 
the core of the galaxy. A spherically symmetrical star form- 
ing population would naturally develop in the centre of the 
galaxy, and one would not need the extra ingredient of a 
short relaxing time-scale as in the merger scenario. 

The central blue cores arise naturally in our models, in 
this way reproducing the colour gradients which are typical 
of BCSs. Within the scenario of our model, if the HDF-2.264 
BCS at z=0.478 is typical of BCSs, then they are z ~ 0.5 
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counterparts of sub-L* spheroids, formed at redshifts some- 
what larger than zgf = 1, and requiring very prolonged 
(> 3 Gyr) cooling flow- fed central stax formation or recur- 
rent central star formation episodes. This could account for 
the scarcity of BCSs in rich clusters due to the removal of 
the gas in the dark halo of the spheroid by the intraclustcr 
medium. Without the gas reservoir in the dark halo, the gas 
inside the galaxy would escape very rapidly, thus preventing 
any long term cooling flow that could feed very prolonged 
central star formation activity. At the same time, any late 
central starburst would be suppressed, since it could not 
be feed either by gas ejected by evolved stars and accumu- 
lated in the galaxy or by old galactic wind ejecta braked in 
the gas-filled dark matter halo and falling back towards the 
galaxy. 

If the progenitors of ~ L* spheroids were formed at 
redshifts larger than zgf = 1, they would give rise to high- 
mass, high-redshift counterparts of the 0.4 < z < 1 BCSs. 
If they arc formed at zgf ~ 5, they could be missed in 
the standard multi-colour surveys designed to select candi- 
dates to LBGs. Although it is possible to find some BCSs in 
the present LBG samples, the colour criteria used in search- 
ing LBGs should be made more flexible (for instance, red- 
der -B450 — Vfjoe colour indices) in order to search for high- 
redshift BCSs which would be the 0.3— a to few Gyr old 
progenitors of the present day ~ L* ellipticals. Moreover, 
according to the chemodynamical model for evolution of 
spheroids, due to the fact that the prolonged central star 
formation characteristic of BCSs, extends at least 2 Gyr, 
and, sometimes, longer than 3 Gyr, the BCS stage in the 
evolution of spheroids lasts somewhat longer than the LBG 
stage (0.2 ^ t ^ 1.5 Gyr) (FT99) and should also be found 
in abundant numbers at high redshifts, provided that we ex- 
plore regions in the colour diagrams wider than those used in 
the selection of LBG candidates, and perform deeper imag- 
ing with higher resolution of the candidates to BCSs. Fi- 
nally, those high-redshift BCSs would be useful probes for 
the epoch of galaxy formation and the cosmological param- 
eters of the Universe. 
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